Go et al. BMC Genomics 2012, 13:627 
http://www.biomedcentral.eom/1 471 -21 64/1 3/627 



^BMC 



Genomics 



RESEARCH ARTICLE Open Access 



2009 pandemic H1N1 influenza virus elicits 
similar clinical course but differential host 
transcriptional response in mouse, macaque, 
and swine infection models 

Jennifer T Go 1 " 1 , Sarah E Belisle 11 , Nicolas Tchitchek 1 , Terrence M Tumpey 4 , Wenjun Ma 3 , Juergen A Richt 3 , 
David Safronetz 5 , Heinz Feldmann 5 and Michael G Katze 1,2 



Abstract 

Background: The 2009 pandemic H1N1 influenza virus emerged in swine and quickly became a major global 
health threat. In mouse, non human primate, and swine infection models, the pH1N1 virus efficiently replicates in 
the lung and induces pro-inflammatory host responses; however, whether similar or different cellular pathways 
were impacted by pH1N1 virus across independent infection models remains to be further defined. To address this 
we have performed a comparative transcriptomic analysis of acute phase responses to a single pH1N1 influenza 
virus, A/California/04/2009 (CA04), in the lung of mice, macaques and swine. 

Results: Despite similarities in the clinical course, we observed differences in inflammatory molecules elicited, and 
the kinetics of their gene expression changes across all three species. We found genes associated with the retinoid 
X receptor (RXR) signaling pathway known to control pro-inflammatory and metabolic processes that were 
differentially regulated during infection in each species, though the heterodimeric RXR partner, pathway associated 
signaling molecules, and gene expression patterns varied among the three species. 

Conclusions: By comparing transcriptional changes in the context of clinical and virological measures, we 
identified differences in the host transcriptional response to pH1N1 virus across independent models of acute 
infection. Antiviral resistance and the emergence of new influenza viruses have placed more focus on developing 
drugs that target the immune system. Underlying overt clinical disease are molecular events that suggest 
therapeutic targets identified in one host may not be appropriate in another. 
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Background 

In 2009 a novel H1N1 influenza virus emerged and rap- 
idly spread worldwide [1]. Clinical disease among 
affected individuals ranged from mild infection to more 
severe pneumonia associated with acute respiratory dis- 
tress syndrome [2,3]. A recent study estimates that over 
284,000 deaths occurred globally within the first year of 
the pandemic [4]. Highly pathogenic influenza virus in- 
fection is characterized by a powerful and potentially 
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destructive immune response [5,6]. Enhanced cytokine 
production has been observed in cynomolgus macaques 
infected with 2009 pandemic H1N1 influenza (pHlNl) 
virus compared to former seasonal H1N1 influenza virus 
[7,8], though pHlNl virus does not precipitate a cyto- 
kine storm' seen with highly pathogenic influenza viruses 
(reviewed in [9]). 

The pHlNl virus has been intensively studied in mice, 
macaques and swine, among other animal models of 
influenza infection. In the absence of prior host adap- 
tation, inoculation of mice with pHlNl virus generally 
results in nonlethal infection that is resolved around day 
8 post-infection. There is low morbidity (<10% total body 
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weight) and moderate bronchiolitis observed in the lung, 
despite efficient viral replication throughout infection 
[10]. Alteration of viral PB1-F2 and NS1 proteins mar- 
ginally impacted viral pathogenesis [11,12] and did not 
affect viral control of innate immune responses during 
pHlNl infection in mice [13]. Cynomolgus macaques 
infected with pHlNl virus have shown diffuse alveolar 
damage, more severe pulmonary lesions, and efficient 
replication of the virus in the lungs compared to seasonal 
H1N1 influenza virus [7,8,14]. Swine infected with 
pHlNl virus develop mild respiratory disease character- 
ized by coughing, sneezing, and acute bronchiolitis, with 
active replication in the lung [15]. Recent studies have 
shown that the Eurasian-origin NA and M genes of 
pHlNl virus are important for transmission among 
swine [16] and can also confer transmission in ferrets 
[17] and guinea pigs [18]. 

Differential gene expression induced by pHlNl virus 
compared to concurrent human or swine influenza 
viruses has been investigated in mice, macaques and 
swine by microarray. In mice, the transcriptional pro- 
gram suggests that regulation of lipid metabolism genes 
plays a protective role, differentiating non-lethal wild- 
type and lethal mouse-adapted pHlNl virus infection 
[10]. Macaques infected with clinical pHlNl isolates 
from Mexico showed increased expression of NFkB 
signaling molecule genes, cytokine and chemokine genes, 
as well as antigen presentation pathway genes compared 
to seasonal H1N1 influenza virus [7]. In a separate study, 
macaques infected with pHlNl virus also showed 
enhanced expression of antiviral and interferon (IFN)- 
regulated genes, such as IFIT2 and ISG15, compared to 
seasonal influenza virus [14]. Transcriptomic analysis of 
swine host responses to pHlNl virus revealed pro- 
nounced inflammatory response gene expression accom- 
panied by increased expression of PPARG-associated 
lipid metabolism genes compared to the 1918-like classical 
swine influenza A/swine/Iowa/ 15/ 1930 (H1N1) virus [15]. 
In contrast to in vivo models, transcriptomic profiling of 
infected type I alveolar epithelial cells showed comparable 
IFN-mediated antiviral and cytokine responses to pHlNl 
and seasonal H1N1 influenza viruses [19]. 

Performing a comparative transcriptomic analysis 
across multiple model systems has revealed conserved 
responses during influenza infection, as demonstrated in 
a recent study by McDermott et al. that used multivari- 
ate modeling approaches to identify similarities in tran- 
scriptional responses to H5N1 virus in the lungs of mice 
and macaques, and human lung epithelial cells [20]. To 
gain further insight into host responses to pHlNl virus 
during acute infection, we examined lung gene expres- 
sion from mice, macaques and swine infected with 2009 
pandemic H1N1 influenza A/California/04/2009 (CA04) 
virus and compared the transcriptional response in each 



host. Our goals were to identify shared or differential 
gene expression patterns across species and to infer 
potential regulators mediating these changes during 
acute infection. Even though CA04 virus elicited a similar 
clinical outcome in each of these species, we found sig- 
nificant differences in the expression of inflammatory 
and lipid metabolism genes, likely impacted by nuclear 
hormone receptor signaling complexes including LXR/ 
RXR that is known to regulate cholesterol homeostasis 
during inflammation. A greater understanding of the 
differences in acute responses from different hosts is 
important because it will aid in the design of tailored 
immunotherapies to influenza virus. 

Results and discussion 

We have previously shown 2009 pandemic H1N1 influ- 
enza A/California/04/2009 (CA04) virus efficiently repli- 
cates in the lung of mice and swine, inducing expression 
of pro-inflammatory genes and causing acute bronchiolitis 
[10,15]. Infection of macaques with CA04 virus is marked 
by productive virus replication in the respiratory tract 
and moderately severe clinical symptoms peaking on day 

6 post-infection (p.i.), with resolution by day 14 p.i. 
( a Safronetz and Feldmann, personal communications). 
Here, we examined mouse, macaque and swine responses 
to CA04 virus by microarray to further investigate tran- 
scriptional changes during acute pHlNl infection. Mice 
were inoculated with 10 6 PFU virus and three animals 
were euthanized on days 1, 3, and 5 p.i. These three time 
points were selected for studying acute phase responses in 
the mouse. Cynomolgus macaques were inoculated with a 
total infectious dose of 7xl0 6 TCID 50 and two animals 
were euthanized on days 1 and 6 p.i. Day 1 macaque lung 
samples were collected for studying early host immune 
responses corresponding to early disease progression and 
day 6 lung samples were chosen for maximum pathology. 
Swine were inoculated intratracheally with 10 6 TCID 50 / 
animal and five animals were euthanized on days 3, 5 and 

7 p.i. Swine lung samples were collected on days 3 and 5 
to examine acute phase responses corresponding to max- 
imal virus shedding and day 7 lung samples were chosen 
for recovery phase. Due to differences in the timing and 
kinetics of acute phase responses in each animal model, 
we examined gene expression changes within each host 
and then compared across species. 

Viral mRNA expression was measured in each sample 
to verify CA04 infection prior to microarray analysis 
(Additional file 1: Figure SI). In general, viral mRNA 
levels correlated with viral titers measured in the lung 
for each species. In mice, viral mRNA expression averaged 
between 3.6 - 4.5 log 10 RQ on day 1 to 5 p.i., correspond- 
ing to average lung virus titers between 5.2 - 6.1 log 10 
PFU/g of tissue on day 1 to 5 p.i. [10]. In macaques, aver- 
age viral mRNA expression was highest on day 1 and 
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decreased by day 6 pi, which was also reflected in the 
virus lung titers that reached around 10 7 and 10 4 TCID 50 
(log 10 )/g on days 1 and 6 p.i., respectively ( a Safronetz and 
Feldmann, personal communications). In swine, viral 
mRNA levels paralleled virus shedding (from days 3 to 5 
p.i.), and the decreased viral mRNA expression observed 
on day 7 p.i. was concomitant with resolution of infection 
in these animals [15]. Quantitation of viral mRNAs in 
infected samples used for microarray allowed us to correl- 
ate host gene expression changes with relative infection 
levels. 

Gene expression analysis of mouse, macaque and swine 
lung infected with CA04 virus 

Gene expression was profiled using species-specific 
commercial oligonucleotide arrays. Each array contained 
a different set of transcripts, with greater redundancy for 
macaque and swine compared to mouse. There are a 
total of 43,379 probes represented on the Mouse Whole 
Genome Gene Expression Microarray, a total of 20,217 
probes represented on the Rhesus Macaque Gene Ex- 
pression Microarray, and a total of 47,813 probes repre- 
sented on the Porcine Gene Expression Microarray VI. 
Due to differences in gene annotation across species, dif- 
ferential gene expression analysis was performed for the 
transcripts associated with 4118 unique genes common 
to all three arrays using Ensembl gene identifiers asso- 
ciated with each array probe. This strategy accounted for 
differences in the number of probes with annotated tran- 
scripts and reduced the likelihood of falsely identifying 
differences in gene expression due to gaps in annotation 
or gene representation. Within this set of genes, there 
was a large representation of genes associated with Cell 
Death, Cancer, and Cellular Growth and Proliferation 
functional categories, as well as genes associated with 
Glucocorticoid Receptor Signaling, IL-12 Signaling and 
Production in Macrophages, and Acute Phase Response 
Signaling canonical pathways. 

We investigated the functional classes of underrepre- 
sented gene sets present on each specie-specific array, as 
the focus on unique genes common to all three arrays 
may introduce bias in our analyses. There was enrich- 
ment of Molecular Mechanisms of Cancer, Axonal 
Guidance Signaling and G-Protein Coupled Receptor 
Signaling canonical pathways, and Gene Expression 
functional annotations related to transcription, organis- 
mal death and abnormal morphology of cells. Data 
integration and interpretation with a cross-species tran- 
scriptomic analysis brings its own challenges and al- 
though complete physical maps have been developed for 
mouse, macaque and swine genomes to support genome 
sequencing and comparative genomics, functional anno- 
tation to date is mostly based on human, mouse, and rat 
literature. As annotation improves, particularly for less 



characterized species such as swine, we will likely be 
able to more fully understand host responses to influ- 
enza virus using microarray and next-generation sequen- 
cing technologies. 

Of the 4118 unique genes common to all three species 
arrays, we identified a total of 696 differentially 
expressed (DE) genes in mice, 771 DE genes in macaque, 
and 611 DE genes in swine that significantly changed in 
response to CA04 virus on at least one day (Students 
t-test P < 0.05, average fold-change > 2). The Venn 
diagram shown in Figure 1 illustrates the overlap of 
CA04 virus-induced DE genes among mouse, macaque 
and swine infection models, and only 53 genes were 
differentially expressed in all three species. Functional 
analysis of each specie DE gene set revealed significant 
enrichment of genes associated with Acute Phase Re- 
sponse Signaling, LXR/RXR Activation, and Athero- 
sclerosis Signaling canonical pathways in infected mice, 
VDR/RXR Activation, Hematopoiesis from Pluripotent 
Stem Cells, and Communication between Innate and 
Adaptive Immune Cells canonical pathways in infected 
macaques, and LXR/RXR Activation, Nitrogen Metab- 
olism, and Antigen Presentation Pathway canonical 
pathways in infected swine (Table 1). The genes asso- 
ciated with each canonical pathway shown in Table 1 
are reported in Additional file 2. 

Inflammatory Response was significantly enriched in all 
three CA04 infection models, with predicted increased ac- 
tivation in macaques and swine (corrected z-scores of 
3.437 and 4.02, respectively) (Additional file 3: Table S2). 

Macaque 
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242 98 207 




Figure 1 Comparison of significant host genes differentially 

regulated during CA04 infection in each species. Venn diagram 

showing the overlap between mouse (red), macaque (green) and 

swine (blue) DE gene sets. 
\ ) 
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Table 1 Canonical Pathways enriched in mouse, macaque and swine CA04 infection models 


Infection model 


Canonical pathway 


B-H p-value 


mol a /mol b 


Mouse 


Acute Phase Response Signaling 


2.34E-10 


28/177 




LXR/RXR Activation 


2.34E-10 


26/136 




Atherosclerosis Signaling 


3.47E-09 


22/129 




Hepatic Fibrosis / Hepatic Stellate Cell Activation 


6.31 E-07 


21/147 




IL-12 Signaling and Production in Macrophages 


1.12E-06 


20/155 




Complement System 


4.07E-06 


10/35 




Type I Diabetes Mellitus Signaling 


5.50E-06 


17/121 




Nitrogen Metabolism 


8.51E-05 


9/120 




Clathrin-mediated Endocytosis Signaling 


1 .30E-04 


20/195 




Dendritic Cell Maturation 


1 .70E-04 


18/185 


Macaque 


VDR/RXR Activation 


3.63E-03 


13/81 




Hematopoiesis from Pluripotent Stem Cells 


5.25E-03 


10/64 




Communication between Innate and Adaptive Immune Cells 


2.75E-02 


12/109 




Type I Diabetes Mellitus Signaling 


3.16E-02 


13/121 




Altered T Cell and B Cell Signaling in Rheumatoid Arthritis 


3.80E-02 


11/92 




Gal 2/1 3 Signaling 


3.80E-02 


13/127 




Neuroprotective Role of TH0P1 in Alzheimer's Disease 


3.80E-02 


7/54 




Hematopoiesis from Multipotent Stem Cells 


3.80E-02 


4/12 




Primary Immunodeficiency Signaling 


3.80E-02 


8/63 




Pyrimidine Metabolism 


3.80E-02 


14/215 


Swine 


LXR/RXR Activation 


1.15E-07 


20/136 




Nitrogen Metabolism 


9.33E-04 


8/120 




Antigen Presentation Pathway 


9.33E-04 


8/40 




Communication between Innate and Adaptive Immune Cells 


9.33E-04 


12/109 




Hepatic Fibrosis / Hepatic Stellate Cell Activation 


9.33E-04 


15/147 




Type I Diabetes Mellitus Signaling 


9.33E-04 


13/121 




Primary Immunodeficiency Signaling 


9.33E-04 


9/63 




Role of Pattern Recognition Receptors in Recognition of Bacteria and Viruses 


9.33E-04 


12/106 




TREM1 Signaling 


9.33E-04 


9/66 




Atherosclerosis Signaling 


1.12E-03 


13/129 



Ingenuity Pathway Analysis was used to determine the top 10 Canonical Pathways. Benjamini-Hochberg (B-H) Multiple Testing Correction p-value was used to 

rank the significance associated for each pathway. 

a Number of molecules differentially expressed in the Canonical Pathway. 

b Total number of molecules in the annotated Canonical Pathway. 



A similar number of DE inflammatory response genes 
characterized CA04 infection in mice, macaques and 
swine (175 DE genes in mice; 166 DE genes in macaques; 
and 160 DE genes in swine), though the kinetics of the 
host response was distinct in each species (Figure 2). In 
general, mice showed an increase in gene expression as in- 
fection progressed, whereas the majority of these genes 
were consistently upregulated or downregulated on days 1 
and 6 pi in CA04 virus-infected macaques. In contrast to 
mice and macaques, swine exhibited enhanced gene ex- 
pression on days 3 and 5 p.i. and the host response 
tapered as the infection resolved. The dissimilarity in the 
kinetics of gene expression across species may be in part 



due to differences in virus replication and timing of acute 
phase responses. Therefore, a more complete kinetic time 
course, including a range of inoculation doses, would be 
necessary to fully appreciate the impact of differences in 
acute phase responses to pHlNl virus infection. 

Differential inflammatory response gene expression 
distinguishes CA04 virus infection in mice, macaques and 
swine 

Among the inflammatory response genes there was rep- 
resentation of Acute Phase Response Signaling, LXR/ 
RXR Activation, and Dendritic Cell Maturation canon- 
ical pathways in infected mice, Role of NFAT in 
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Mouse Macaque Swine 

-0.6 0 +0.6 

Figure 2 Heatmap of inflammatory response genes across all 
three species infected with CA04 virus. Average log 10 (ratio) 
expression of 319 DE inflammatory response genes induced by 
CA04 virus. In mice and swine, infected lung gene expression is 
referenced to specie-matched mock at each time point. In 
macaques, infected lung gene expression is relative to an uninfected 
lung reference pool at each time point. Red indicates expression 
was increased relative to the control reference and green indicates 
that expression was decreased relative to the control reference. 
Saturation is 4-fold. 



Regulation of the Immune Response, Role of Macro- 
phages, Fibroblasts and Endothelial Cells in Rheumatoid 
Arthritis, and G-Protein Coupled Receptor Signaling ca- 
nonical pathways in infected macaques, and enrichment 
of inflammatory response genes associated with Gluco- 
corticoid Receptor Signaling, Dendritic Cell Maturation, 
and LXR/RXR Activation canonical pathways in infected 
swine. When we examined the 53 genes common to 
mice, macaques and swine, we noted the majority of 
genes were upregulated in response to infection in mice 
on days 3 and 5 p.i., in macaques on days 1 and 6 p.L, 
and in swine on days 3 and 5 p.i. (Figure 3). Inflamma- 
tory response genes (highlighted yellow) included inter- 
feron (IFN) signaling molecule, STAT1, IFN-regulated 
antigen presentation and immunoproteasome compo- 
nents, TAP1, PSMB8 and PSMB9, Toll-like receptor 4 
coreceptor, CD 14, and antiviral effectors, IFIH1 and 
IFIT2, previously shown to play essential roles in innate 
immune response to influenza infection. We also found 
several leukocyte-specific genes commonly differentially 
expressed during CA04 infection. Members of the 



immunoglobulin receptor superfamily upregulated 
across species included CD72, CD274 (also known as 
programmed cell death ligand 1), and CD180, a TLR 
homologue expressed on B cells, macrophages and den- 
dritic cells. Macrophage-restricted receptor, SIGLEC1, 
neutrophil factors, PLUNC and NCF4, and T lympho- 
cyte coreceptor, CD8A, were also differentially expressed 
during CA04 infection in each species (Figure 3, 
Additional file 4: Table S3). 

Several models show an association between enhanced 
immune cell infiltrate and severe lung immunopathol- 
ogy. Excessive macrophage and neutrophils are observed 
in the lung of mice following H5N1 and 1918 infection 
[21], and in pregnant animals infected with pHlNl virus 
[22]. In mice, macaques, and swine infected with CA04 
virus, leukocyte and lymphocyte responses were evident 
in all three species based on immune cell-specific gene 
expression changes detected in the lung, though the 
gene expression patterns varied across species. For ex- 
ample, there was greater upregulation of macrophage 
factor, SIGLEC1, observed in mice (days 3 and 5 p.i.) 
and macaques (days 1 and 6 p.i.) compared to swine on 
day 3 p.i. Enhanced expression of leukocyte associated 
immunoglobulin-like receptor 1, LAIR1, was observed in 
macaques and swine infected with CA04 virus, as com- 
pared to infected mice. In contrast to mice and maca- 
ques, swine exhibited strong expression of neutrophil 
factor, PLUNC, which we noted was downregulated dur- 
ing infection in the other species (Figure 3). These 
results reflect a rapid inflammatory shift in the lungs of 
mice, macaques, and swine during CA04 infection that 
involves different immune cell responses, but unlike 
H5N1 and 1918 viruses, these immune cell responses do 
not cause immunopathology. A cross-species transcrip- 
tomic comparison of the host response to H5N1 and 
1918 viruses would be necessary to further explore the 
potential role of specific immune cell responses to influ- 
enza pathogenesis. 

CA04 infection impacts glucocorticoid receptor signaling 
differently across species 

We further sought to evaluate CA04 virus-induced host 
responses that were unique to each species by investigat- 
ing non-overlapping gene sets shown in Figure 1. The 
most significant canonical pathways represented in each 
of these three gene sets are summarized in Additional 
file 5: Figure S2. The 207 DE gene set unique to swine 
was enriched for genes associated with Role of JAK1 and 
JAK3 in yc Cytokine Signaling canonical pathway that 
included JAK3 and STAT5A genes. Differential regula- 
tion of JAK3/STAT5 signaling in CA04 virus-infected 
swine may impact the homeostasis and activation of per- 
ipheral T lymphocytes. For example, STAT5A/B defects 
results in enhanced apoptosis of T lymphocytes in mice 
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Figure 3 Heatmap of 53 genes commonly differentially regulated across all three species infected with CA04 virus. Average log 10 (ratio) 
expression of 53 DE genes induced by CA04 virus. In mice and swine, infected lung gene expression is referenced to specie-matched mock at 
each time point. In macaques, infected lung gene expression is relative to an uninfected lung reference pool at each time point. Red indicates 
expression was increased relative to the control reference and green indicates that expression was decreased relative to the control reference. 
Saturation is 4-fold. A total of 22 DE genes associated with inflammatory response are highlighted on the right; yellow denotes positive 
association. 



[23]. Differential expression of JAK3 and STAT5A genes 
in swine during CA04 infection may also suggest a bal- 
ance between inflammatory host defenses and gluco- 
corticoid receptor (GR)-mediated cellular growth and 
survival. Closer inspection of each species inflammatory 
response DE gene set showed noteworthy representation 
of GR Signaling related genes, such as A2M, FGG, 
HSPA5, IL4, and MAPK13 genes in mice, FKBP4, JAK1, 
NCOA2 (also known as GRIP-1), and NRAS genes in 
macaques, and HSPA4 (also known as HSP70) in 
addition to JAK3, and STAT5A genes in swine. Alpha-2- 
macroglobulin (A2M) identified by mass spectrometry in 
human saliva was found to exhibit antiviral activity 
against pHlNl virus [24]. 



Since GR signaling appears to play a role during CA04 
infection in mice, macaques, and swine, we used a tran- 
scription factor (TF) analysis strategy that combined 
transcriptomic and DNA sequence information to inves- 
tigate GR and other potential regulators controlling host 
responses during CA04 infection (described in Meth- 
ods). As summarized in Table 2, we found predicted in- 
hibition of hepatocyte nuclear factor lalpha, HNF1A, in 
infected mice on all three days, consistent with our pre- 
vious report [10], and predicted activation of pro- 
inflammatory regulators, STAT1, NFkB, and IRF7 in at 
least two species. Examination of downstream gene tar- 
gets of GR determined using IPA Upstream Regulator 
Analysis showed differences in expression across species 
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Table 2 Upstream Regulator analysis of mouse, macaque and swine CA04 infection models 



Upstream regulator 


Molecule type 


p-value of overlap 


Regulation z-score (dpi) 


Predicted regulator status 


Mouse 










STAT3 


transcription regulator 


1.68E-17 






STAT1 


transcription regulator 


2.80E-17 


2.747 (3); 2.124 (5) 


Activated 


PPARA 


ligand-dependent nuclear receptor 


1.12E-15 






NFkB (complex)* 


complex 


2.37E-13 


2.463 (3); 2.787 (5) 


Activated 


JUN 


transcription regulator 


6.74E-13 


2.072 (5) 


Activated 


YY1 


transcription regulator 


4.65 E- 12 






SMAD3 


transcription regulator 


5.13E-12 






IRF7 


transcription regulator 


3.52E-1 1 


2.994 (3); 2.532 (5) 


Activated 


IRF1 


transcription regulator 


4.14E-11 


2.125 (3) 


Activated 


HNF1 A* 


transcription regulator 


5.01 E-11 


-3.166(1); -3.743 (3); -3.035 (5) 


Inhibited 


Macaque 










IRF1 


transcription regulator 


2.16E-10 


2.547 (1) 


Activated 


SPI1 


transcription regulator 


4.74E-08 


2.490 (1) 


Activated 


TRIM24 


transcription regulator 


2.43E-07 


-3.395 (1,6) 


Inhibited 


IRF7 


transcription regulator 


4.98E-07 


4.284(1,6) 


Activated 


STAT3 


transcription regulator 


4.31E-06 


2.096 (1,6) 


Activated 


STAT1 


transcription regulator 


1.15E-05 


3.232 (1); 3.57 (6) 


Activated 


JUN 


transcription regulator 


1.57E-05 






IRF8 


transcription regulator 


7.08E-05 






VitaminD3-VDR-RXR 


complex 


7.75E-05 






ELF1 


transcription regulator 


9.71 E-05 






Swine 










STAT1* 


transcription regulator 


6.54E-14 


2.823 (3); 3.157 (5) 


Activated 


TRIM24 


transcription regulator 


2.59E-1 1 






IRF1* 


transcription regulator 


7.54E-1 1 






PPARA 


ligand-dependent nuclear receptor 


2.34E-09 


-2.001 (3) 


Inhibited 


STAT3 


transcription regulator 


2.35E-09 






NFkB complex 


complex 


2.48E-09 


3.01 (5) 


Activated 


IRF2* 


transcription regulator 


7.64E-09 






IRF7 


transcription regulator 


1 .52E-08 






JUN 


transcription regulator 


4.38E-08 






NFKB1 


transcription regulator 


2.32E-07 


2.136 (7) 


Activated 



Ingenuity Pathway Analysis was used to determine the top 10 Upstream Regulators. The p-value of overlap was used to rank the significance associated for each 
Upstream Regulator. The p-value indicates the significance of the overlap between the genes targeted by the upstream regulator in the IPKB database and the 
experimental dataset. Genes with an apteryx signify regulators that were identified as significantly enriched (P < 0.0001) by PSCAN. Z-scores for predicted 
upstream regulators (|z|>2) in each species and at each time point are shown. Z > 2 predicts activation of the upstream regulator. Z < -2 predicts inhibition of the 
upstream regulator. 



with more pronounced transcriptional changes observed 
in mice on days 3 and 5 p.i. and on the early time points 
in macaques and swine (Figure 4). Through a DNA 
sequence-based approach, we assessed the representa- 
tion of DNA-binding motifs in the promoter region of 
DE genes targeted by HNF1A and STAT1 using DNA 
sequences from each species genome and regulator 
DNA-binding preferences, modeled as affinity Position 
Weight Matrices (PWM) (Additional file 6: Figure S3). 



We have identified enriched TFs and targets from our 
dataset; however, the analyses reported here coordinately 
investigate the expression and sequence-specific TF 
DNA binding sites of the DE genes, and it does not in- 
vestigate the combinatorial effects of multiple factors 
binding to the genomic regions or molecular determi- 
nants of transcriptional responses such as receptor 
ligation or histone acetylation that can collectively con- 
tribute to distinct transcriptional profiles [25]. In 
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Figure 4 Glucocorticoid receptor target genes are differentially 
expressed in response to CA04 virus. Average log 10 (ratio) 
expression of downstream GR target genes induced by CA04 virus. 
In mice and swine, infected lung gene expression is referenced to 
specie-matched mock at each time point. In macaques, infected 
lung gene expression is relative to an uninfected lung reference 
pool at each time point. Red indicates expression was increased 
relative to the control reference and green indicates that expression 
was decreased relative to the control reference. Saturation is 4-fold. 



addition, restriction of the search space to -450 to +50 
nucleotides relative to the transcription start site (TSS) 
can preclude the identification of more distal TFs effect- 
ing expression, which may explain why GR is not identi- 
fied as a predicted regulator (Table 2), although many 
target genes are differentially expressed in the data set 
(Figure 4). Future investigations of GR-mediated host 
responses to CA04 virus would need to take into ac- 
count these considerations as well as the physical inter- 
actions between GR and STAT5, for example, and 
chromatin modifications known to occur [26] to provide 
a more comprehensive view of transcriptional regulation 
during CA04 infection across different species. 

CA04 virus significantly alters expression of genes 
involved in cholesterol homeostasis in mice and swine, 
and vitamin D receptor genes in macaques 

Vitamin D receptor (VDR) and liver X receptor (LXR) 
are ligand- activated transcription factors of the nuclear 
receptor superfamily that heterodimerize with retinoid X 
receptor (RXR) to modulate an array of immune and 
metabolic programs (reviewed in [27,28]). Perturbation 
of RXR-mediated signaling pathways have shown unique 
roles for VDR in macrophage responses to Mycobacter- 
ium tuberculosis [29], and LXR in macrophage responses 
to Listeria monocytogenes [30]. In mice and swine 
infected with CA04 virus, we found LXR/RXR Activa- 
tion was among the most significant canonical pathways 
differentially regulated during acute infection, while 
VDR/RXR Activation was the most significant canonical 
pathway differentially regulated in CA04 virus-infected 
macaques (Table 1). As shown in Figure 5, gene expres- 
sion profiles of LXR/RXR and VDR/RXR associated 
pathway molecules, as well as genes enriched for Com- 
munication between Innate and Adaptive Immune Cells 
and Acute Phase Response Signaling canonical pathways 
strongly emphasize the distinct nature of the host re- 
sponse to CA04 virus in these three animal models. 

Macrophages have been shown to increase in mouse 
lungs following influenza infection and these cells are 
specific targets of viral infection [21,22]. In response to 
CA04 virus, we observed a number of genes related to 
macrophages and with roles in lipid metabolism, such as 
cholesterol efflux, when examining LXR/RXR Activation 
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Figure 5 Immune genes exhibit unique expression profiles in mice, macaques, and swine infected with CA04 virus. Average log 10 (ratio) 
expression of DE genes induced by CA04 virus. In mice and swine, infected lung gene expression is referenced to specie-matched mock at each 
time point. In macaques, infected lung gene expression is relative to an uninfected lung reference pool at each time point. Red indicates 
expression was increased relative to the control reference and green indicates that expression was decreased relative to the control reference. 
Saturation is 4-fold. There are a total of 16 DE genes represented for Communication between Innate and Adaptive Immune Cells, 28 DE genes 
represented for Acute Phase Response Signaling, 31 DE genes represented for LXR/RXR Activation, and 12 DE genes represented for VDR/RXR 
Activation. The association of each DE gene within a given species for each of the represented Canonical Pathways is shown on the left; blue 
denotes differential expression and gray denotes no differential expression. The following abbreviations are used: Mm, Mus musculus (mouse); 
Mf, Macaco fascicularis (macaque); and Ss, Sus scrota (swine). 



in mice and swine. In mice, there was marked decreased 
expression of apolipoprotein genes, including APOE and 
APOA1, and moderate downregulated expression of 
macrophage-specific genes involved in recognition of 
oxidized lipids, such as macrophage scavenger receptor 
genes MSR1 (SR-A1) and CD36 (Figure 6, top panel; 
Additional file 7: Table S4). The ABCA1 gene, known 
for its role in cholesterol transport, and the LPL gene, 
involved in lipid catabolism, were also downregulated in 
mice. In contrast to mice, swine showed strong up- 
regulation of these genes associated with cholesterol ef- 
flux on day 3 p.i. Previous studies have shown a 
decrease in ABCA1 mRNA expression in LXR-activated 
primary macrophages via an IRF3-dependent mechanism 



following influenza infection [31]. Here, we observe that 
ABCA1 and other genes associated with cholesterol ef- 
flux are significantly downregulated in response to CA04 
infection in mice, but upregulated in CA04 virus - 
infected swine lung, suggesting that mice may be more 
efficient at suppressing cholesterol efflux, and in turn, 
possibly preventing the formation of lipid-laden macro- 
phages (also known as foam cells). Foam cells are the 
product of inflammatory responses ranging from athero- 
sclerosis (reviewed in [32]) to infection by Chlamydia 
pneumonia [33] and Toxoplasma gondii [34]. While the 
link between influenza-induced innate immunity and 
cholesterol metabolism is poorly defined, our findings 
suggest a role for LXR/RXR signaling during CA04 
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Figure 6 (See legend on next page.) 
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(See figure on previous page.) 

Figure 6 CA04 virus induces different RXR signaling-mediated responses in each species. A) Average log 10 (ratio) expression of 31 LXR/RXR 
DE genes induced by CA04 virus. In mice and swine, infected lung gene expression is referenced to specie-matched mock at each time point. 
In macaques, infected lung gene expression is relative to an uninfected lung reference pool at each time point. Red indicates expression was 
increased relative to the control reference and green indicates that expression was decreased relative to the control reference. Saturation is 
4-fold. B) Average log 10 (ratio) expression of 12 VDR/RXR DE genes induced by CA04 virus. In mice and swine, infected lung gene expression is 
referenced to specie-matched mock at each time point. In macaques, infected lung gene expression is relative to an uninfected lung reference 
pool at each time point. Red indicates expression was increased relative to the control reference and green indicates that expression was 
decreased relative to the control reference. Saturation is 4-fold. The association of each DE gene within a given species is shown on the left; blue 
denotes differential expression and gray denotes no differential expression. The following abbreviations are used: Mm, Mus musculus (mouse); 
Mf, Macaca fascicularis (macaque); and Ss, Sus scrota (swine). 

V J 



infection and indicate that this cellular pathway is 
impacted differently in two distinct hosts. 

In human macrophages and respiratory epithelial cells, 
vitamin D plays a role in innate immune responses by 
increasing TLR coreceptor, CD 14, and stimulating ex- 
pression of antimicrobial peptides, such as cathelicidin 
during respiratory infection [35,36]. Vitamin D has also 
been investigated in epidemiological studies of reduced 
risk of influenza infection in different human patient 
cohorts [37]. In contrast to mice and swine, IL2 gene ex- 
pression was more strongly induced in infected maca- 
ques on both days. This may indicate macaques mount 
enhanced Thl cell responses during CA04 virus infec- 
tion compared to mice and swine (Figure 6, lower 
panel). As shown in Figure 6, we also found increased, 
expression of VDR/RXR Activation associated genes 
such as CYP24A1, encoding a hydroxylase involved in 
vitamin D catabolism, RXRG, and HR, a transcriptional 
corepressor of vitamin D receptor. Intriguingly, anti- 
microbial peptide, CAMP (also known as LL-37), was 
found to be downregulated in mouse and macaques 
infected with CA04 virus, and the decreased expression 
of CAMP in CA04 virus -infected macaques may be 
explained in part by the increased expression of APOA1 
(Figure 6, top panel; Additional file 7: Table S4), which 
has been shown to bind and inhibit CAMP [38]. While 
there are several genes associated with VDR/RXR Acti- 
vation that are differentially expressed in macaques 
infected with CA04 virus, vitamin D control of innate 
immune responses does not appear to play a major role 
during infection. 

Conclusions 

We have shown differences in transcriptional responses 
to a single 2009 pandemic H1N1 influenza virus in three 
independent animal models. Our goal in performing a 
cross-species transcriptomic comparison was to identify 
shared and differential gene expression patterns to better 
evaluate the acute phase in different hosts. We found 
significant differences in expression of inflammatory 
response and lipid metabolism genes, which were likely 
impacted by glucocorticoid receptor and retinoid X 



receptor signaling complexes, such as LXR/RXR that 
is known to regulate cholesterol homeostasis during 
inflammation. The timing and magnitude of the host 
response is critical in determining disease outcome 
[39] and the success of immunomodulatory therapy 
(reviewed in [9]). For example, corticosteroid therapy 
in patients infected with pHlNl virus was found to 
increase the risk of developing severe disease [40]. 
Metagenomic studies using clinical samples will be 
necessary to further understand the host response to 
pHlNl virus and other etiological agents within the 
human population, as recently explored by Greninger 
et al. using a pan-viral microarray and deep sequencing 
to characterize pHlNl virus from human nasopharyn- 
geal aspirates [41]. Despite similar clinical outcomes, 
differences in the host transcriptional response could 
suggest that therapeutic targets identified in one host 
may not be relevant in another. 

The gene expression differences elicited by CA04 virus 
within these three independent animal models demon- 
strate the disparate nature of the host response. There 
are particular advantages and disadvantages in modeling 
disease and immune responses to influenza infection in 
mice, macaques and swine. Mice, though largely resistant 
to infection with most human virus isolates, provide a 
tractable system for genetic manipulation to study key 
genes in the regulation of host responses to influenza 
virus. Non human primates are thought to more closely 
model the human response to influenza virus given their 
genetic and physiological similarities, while swine serve 
as a natural reservoir for influenza A viruses and have 
been linked to the emergence of some of the most not- 
able influenza pandemics in recent history, including the 
H1N1 pandemic in 2009. There has been interest in de- 
termining the extent to which influenza viruses are able 
to cause disease in swine, particularly with regard to 
1918 pandemic influenza virus that does not cause severe 
disease in swine like it does in mice and macaques [42], 
This makes swine a unique model for future systems 
biology analyses that may help to uncover host responses 
contributing toward the emergence or maintenance of 
novel influenza viruses with pandemic potential. 
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Methods 

Virus 

2009 pandemic H1N1 influenza virus A/California/04/ 
2009 (CA04) was isolated from a nasal swab of a < 18 y.o. 
boy from San Diego, California [43]. 

Animal models 

Female BALB/c mice (Mus rnusculus), 6 to 8-week-old, 
were intranasally inoculated with 10 6 plaque-forming 
units of CA04 virus in 50 ul (n = 9) or inoculated with 
50 ul of phosphate-buffered saline (control; n = 8). 
Whole lungs were harvested from infected animals at 
days 1, 3 and 5 post-inoculation (n = 3 per time point) 
and from time-matched control animals (n = 3 on days 
1 and 3 and n = 2 on day 5) for extraction of total RNA 
as previously described [10]. Crossbred pigs (Sus Scrofa), 
4-week-old, were inoculated intratracheally with either 
10 6 TCID 50 /pig egg-derived CA04 virus (n = 15) or 
mock inoculated with non-infectious cell culture super- 
natant (control; n = 15) as described elsewhere [44]. 
Animals were euthanized on 3, 5, and 7 dpi (n = 5 per 
time point). Cynomolgus macaques (Macaca fascicu- 
laris), 4 to 15 y.o., weighing 3.0-8.7 kg, were infected 
with CA04 virus (n = 6) under anesthesia through a 
combination of intratracheal (4 ml), intranasal (0.5 ml 
per nostril), conjunctival (0.5 ml per eyelid) and oral (1 ml) 
routes with a suspension containing lxlO 6 TCID 50 /ml 
(total infectious dose was 7xl0 6 TCID 50 ) as described 
elsewhere [7]. Animals were euthanized on 1 and 6 dpi 
(n = 2 per time point). 

Animal ethics 

Mouse infection experiments were completed at the CDC 
under the guidance of the CDC s Institutional Animal Care 
and Use Committee in an Association for Assessment and 
Accreditation of Laboratory Animal Care International- 
(AAAL AC) -accredited animal facility. Swine infection 
experiments were completed at the Central States Research 
Center (CSRC), Inc BSL-3 facility (Oakland, NE) in com- 
pliance with the CSRC s Institutional Animal Care and Use 
Committee studies. Macaque infection experiments were 
approved by the RML Institutional Animal Care and Use 
Committee (IACUC), and performed following the guide- 
lines of AAALAC by certified staff in an AALAC approved 
facility. 

Real-time qRT-PCR 

Real-time PCR was performed using a Custom TaqMan 
Gene Expression Assay (Applied Biosystems) designed 
for CA04 HA sequence (forward primer: AGCTCAG 
TGTCATCATTTGAAAGGT; reverse primer: GGA 
CATGCTGCCGTTACAC; reporter: TTGGGCCATG 
AACTTG). cDNAs were generated using a QuantiTect 
reverse transcription kit (Qiagen). Samples from 



individual animals were run in quadruplicate. rRNA 
(18S) was used to normalize quantification of the target 
and quantification of normalized target was performed 
using the 2" AACt calculation [45]. HA expression was quan- 
tified relative to expression of an endogenous control for 
each specie sample set that did not change with infection 
and expression in an uninfected lung sample. The following 
TaqMan Gene Expression Assays (Applied Biosystems) 
were used: Mfapla (Assay ID Mm00849648_gH) served as 
the endogenous control for mouse samples, B2M (Assay ID 
Rh02847368_ml) served as the endogenous control for ma- 
caque samples, and RPS6 (Assay ID Ss03374061_gl) served 
as the endogenous control for swine samples. Average 
loglORQ expression is shown for each species at each time 
point ± standard deviation. 

Microarray hybridization 

Total RNA isolated from lung tissue from individual 
animals on each day of euthanasia was used for oligo- 
nucleotide array experiments. For swine and mice, 
RNA isolated from mock-infected animals at each time 
point served as an uninfected reference. For cynomolgus 
macaques, a pool of RNA from the lungs of eight unin- 
fected animals matched for age and sex was used as the 
uninfected reference. NanoDrop ND-1000 and Agilent 
2100 Bioanalyzer instrumentation was used to determine 
the concentration and quality of all RNA samples. A total 
of 1500 ng Cy3-labeled probe was used for microarray 
slide hybridizations, thereby normalizing for the input 
RNA amount. Mouse, macaque and swine samples were 
measured with 4x44K commercial arrays from Agilent 
Technologies designed for each species, Mouse Whole 
Genome Gene Expression Microarray (G4122F), Rhesus 
Macaque Gene Expression Microarray (G2519F; Design 
ID: V2: 026806), and Porcine Gene Expression Micro- 
array VI (G2519F; Design ID: VI: 020109). 

Data normalization and accessibility 

The background corrected data from the Feature Ex- 
traction output were normalized across replicates 
within each species using central tendency nor- 
malization (75% percentile, target 1000) within Gene- 
data (Analyst 7.0). All primary microarray data have 
been deposited in NCBI's Gene Expression Omnibus 
(GEO) under GEO Series (GSE) accession number 
GSE40092. The primary microarray data are also avail- 
able at the University of Washington's Public Micro- 
array Data Download site (http://expression.microslu. 
washington.edu). 

Statistical analysis of microarray data 

For mice and swine, Students t-test was performed on 
background corrected, normalized log-intensity data 
comparing CA04 virus-infected lung gene expression to 
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time- and species -matched mock-infected lung gene 
expression at each time point (Unadjusted P- value < 0.05). 
For cynomologus macaques, Students t-test was per- 
formed on background corrected, normalized log- 
intensity data comparing CA04 virus-infected lung gene 
expression to an uninfected lung reference pool (n = 8) 
matched for age and sex at each time point (Unadjusted 
P- value < 0.05). Differentially expressed genes were then 
filtered to include only genes that changed at least two- 
fold compared to mock on at least one day within each 
species. Statistical comparison using Students t-test with 
the Benjamini-Hochberg multiple testing correction [46] 
(Adjusted P-value < 0.05) and no fold-change parameter 
resulted in very few macaque genes that were statistically 
significant, though there was considerable overlap in the 
DE gene sets identified from the two tests for each 
species. 

Functional analysis of differential gene expression data 

Functional analysis of DE genes was performed using 
Ingenuity Pathways Analysis (Ingenuity Systems), which 
analyzes the experimental dataset in the context of 
known biological response and regulatory networks in 
the Ingenuity Pathways Knowledge Base (IPKB). Ensembl 
human gene annotations were used for functional ana- 
lyses of macaque and swine gene sets (Additional file 8: 
Table S5). The right- tailed Fisher s Exact test was used to 
determine the statistical significance of each biological 
function assigned to the gene expression data, and the 
Benjamini-Hochberg (B-H) Multiple Testing Correction 
was applied to p-values to reduce the likelihood that 
statistical associations were due to random chance. 

Transcription factor enrichment analysis 

Upstream Regulator Analysis in IPA incorporates expres- 
sion of downstream target genes from the experimental 
dataset and compiled knowledge of reported relation- 
ships between regulators and their known target genes 
within IPKB. This analytical tool was used to predict 
upstream regulators and infer their activation state by 
calculating a z-score that determines whether gene expres- 
sion changes for known targets of each regulator (z > 2, 
regulator predicted as "activated" and z < -2, regulator 
predicted as "inhibited". Transcription factor binding 
motif enrichment was performed with PSCAN [47] using 
Position Weight Matrices (PWM) for human and mouse 
species obtained from JASPAR CORE database [48] and 
the promoter of each target gene defined from -450 to 
+50 nucleotides relative to the TSS. PSCAN computes a 
z-test P-value for each regulator, which is an assessment 
of whether there is significant representation (P < 0.001) 
of the regulator DNA-binding motif in promoters of the 
queried genes. 



Transcription factor DNA-binding promoter analysis 

Transcription factor (TF) DNA-binding analysis was per- 
formed using STAT1 (human MA0137.2) and HNF1A 
(mouse MA153.1) PWMs obtained from the JASPAR 
CORE database [48] and STAT1 (mouse M00224) and 
HNF1A (human M00206) PWMs obtained from the 
TRANSFAC database 7.0 [49]. Promoters sequences of 
genes analyzed were defined from -450 to +50 nucleo- 
tides relative to the TSS and were retrieved from the 
Ensembl database (release 66). Macaca mulatta 
sequences were used instead of Macaca fascicularis 
sequences for the TF DNA-binding analysis. P-values 
associated for each target gene were calculated using the 
TFM-Pvalue analytical tool described in [50] and distribu- 
tion of nucleotides amongst the promoter sequences have 
been taken into account for background correction. Sig- 
nificant target genes have been identified as having a 
P-value < 10" 4 

Endnotes 

a Safronetz and Feldmann, personal communications 
Laboratory of Virology, Division of Intramural Research, 
National Institute of Allergy and Infectious Diseases, 
National Institutes of Health, Rocky Mountain Labora- 
tories, Hamilton, MT 59840, USA. 
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Additional file 1: Figure SI. Relative influenza HA gene expression in 
the lung of infected mice, macaques and swine. CA04 virus HA mRNA 
was quantified in each lung sample relative to gene expression of an 
endogenous control that did not change with infection (Mfapla for 
mouse; B2M for macaque; Ss03374061_g1 for swine) and expression in 
an infected lung sample that did not have detectable CA04 HA 
expression. Average loglORQ expression ± SD is shown for each species 
at each time point (mice, n =3 per time point; macaques, n = 2 per time 
point; swine, n = 3 per time point). 

Additional file 2: Table SI. Summary of the top 10 Canonical Pathways 
enriched in mouse, macaque and swine CA04 infection models including 
the differentially expressed genes associated with each pathway. 

Additional file 3: Table S2. Functional analysis of DE genes from mice, 
macaques and swine infected with CA04 virus. 

Additional file 4: Table S3. Gene expression of 53 DE genes commonly 
differentially expressed in mouse, macaque and swine lung infected with 
CA04 virus. 

Additional file 5: Figure S2. Pathway enrichment in mouse, macaque 
and swine infected lung unique to each species. Ingenuity Pathway 
Analysis was used to determine the top 5 Canonical Pathways. Fisher's 
Exact test p-value was used to rank the significance associated for each 
Canonical Pathway. 

Additional file 6: Figure S3. Kinetics of STAT1 and HNF1A target gene 
expression in mouse, macaque and swine infected with CA04 virus. A) 
Transcription factor DNA-binding analysis used matrices obtained from 
JASPAR CORE and TRANSFAC databases. Human matrices were applied 
to Macaca mulatta and Sus Scrofa genome sequence scans. For STAT1, 
human MA0137.2 matrix and mouse M00224 matrix were used. For 
HNF1A, human M00206 matrix and mouse MA153.1 matrix were used. B) 
Average fold change gene expression compared to mock for STAT1 and 
HNF1A target genes identified by genome scans with matrices in part A. 
DE genes are shaded gray and target genes of STAT1 (left panel) and 
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HNF1A (right panel) are shaded purple. The STAT1 DE gene is 
highlighted in red. C) Average log 10 (ratio) expression of STAT 1 (left panel) 
and HNF1A (right panel) target genes identified using IPA Upstream 
Regulator Analysis. In mice and swine, infected lung gene expression is 
referenced to specie-matched mock at each time point. In macaques, 
infected lung gene expression is relative to an uninfected lung reference 
pool at each time point. Red indicates expression was increased relative 
to the control reference and green indicates that expression was 
decreased relative to the control reference. Saturation is 4-fold. 

Additional file 7: Table S4. Gene expression of LXR/RXR and 
VDR/RXR genes in mouse, macaque and swine lung infected with CA04 
virus. 

Additional file 8: Table S5. Summary of CA04 induced DE gene sets 
for each species. 
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